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In the present work, we study the quantum cosmology description of a Friedmann- 
Robertson- Walker model in the presence of a stiff matter perfect fluid and a negative 
cosmological constant. We work in the Schutz's variational formalism and the spatial 
sections have constant negative curvature. We quantize the model and obtain the appro- 
priate Wheeler-DeWitt equation. In this model the states are bounded therefore we com- 
pute the discrete energy spectrum and the corresponding eigcnfunctions. In the present 
work, we consider only the negative eigenvalues and their corresponding eigcnfunctions. 
This choice implies that the energy density of the perfect fluid is negative. A stiff mat- 
ter perfect fluid with this property produces a model with a bouncing solution, at the 
classical level, free from an initial singularity. After that, we use the eigenfunctions in 
order to construct wave packets and evaluate the time-dependent expectation value of 
the scale factor. We find that it oscillates between maximum and minimum values. Since 
the expectation value of the scale factor never vanishes, we confirm that this model is 
free from an initial singularity, also, at the quantum level. 
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2 G. Oliveira-Neto et al. 

1. Introduction 

The great importance of cosmological models where the matter content is rep- 
resented by a stiff matter perfect fluid was recognized since its introduction by 
ZeldovichP^ This perfect fluid has an equation of state of the form, p — aw, with 
a = 1, where w and p are, respectively, the fluid energy density and pressure. It 
can also be described by a massless free scalar field. In order to understand better 
the importance of this perfect fluid for cosmology, one has to compute its energy 
density. In the temporal gauge (N(t) = 1), this quantity is proportional to l/a(i) 6 , 
where N(t) is the lapse function and a{t) is the scale factor. On the other hand, 
in the same gauge, the energy density of a radiative perfect fluid is proportional to 
l/a(t) 4 . This result indicates that there may have existed a phase earlier than that 
of radiation, in our Universe, which was dominated by stiff matter. Due to that im- 
portance, many physicists have started to consider the implications of the presence 
of a stiff matter perfect fluid in FRW cosmological models. The first important im- 
plication of the presence of stiff matter in FRW cosmological models is in the relic 
abundance of particle species Produced after the 'Big Bang' due to the expansion 
and cooling of our Universe. I 3 | 4 | 5 | 6 | 7 | 8 | The 

presence of stiff matter in FRW cos- 
mological models may also help explaining the baryon asymmetry and the density 
perturbations of the right amplitude for the large scale structure formation in our 
Universe r It may also pla y an important role in the spectrum of relic gravity 
waves created during inflation!^ Since there may have existed a phase earlier than 
that of radiation which was dominated by stiff matter some ph ysicists considered 
quantum cosmological models with this kind of matter JHHSIlH 

In the present work, we study the quantum cosmology description of a 
Friedmann-Robertson- Walker (FRW) model in the presence of a stiff matter per- 
fect fluid and a negative cosmological constant. The model has constant negati vely 
curved spatial sections. We work in the Schutz's variational formalism! ^^ We 
quantize the model and obtain the appropriate Wheeler-DeWitt equation. In this 
model the states are bounded therefore we compute the discrete energy spectrum 
and the corresponding eigenfunctions. In the present work, we consider only the neg- 
ative eigenvalues and their corresponding eigenfunctions. This choice implies that 
the energy density of the perfect fluid is negative. A stiff m atter per fect fluid with 
this property has already been considered in the literature! 1 ^ | 1 8 | 1 9] p roc | uces a 
model with a bouncing solution, at the classical level, free from an initial singular- 
ity. After that, we use the eigenfunctions in order to construct wave packets and 
evaluate the time-dependent expectation value of the scale factor. We find that it 
oscillates between maximum and minimum values. Since the expectation value of 
the scale factor never vanishes, we confirm that this model is free from an initial 
singularity, also, at the quantum level. 

The presence of a negative cosmological constant in the present model implies 
that the universe described by it has a maximum size, in other words it is bounded. 
Taking into account the current cosmological observations, a negative cosmological 
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An Early Universe Model with Stiff Matter and a Cosmological Constant 3 

constant will not be able to describe the present accelerated expansion of our Uni- 
verse. It is not our intention to describe the present state of our Universe with this 
model. On the other hand, it is our intention to describe a 'possible' state of our 
primordial Universe. One important theory which is a strong candidate to desc ribe 
the unification of all known physical interactions is superstring theoryPQ12U Due 
to that, many physicists believe that superstring theory will correctly describe the 
quantum gravity effects that took place at the beginning of our Universe. There is 
an important conjecture which tells that Type IIB string theory on (AdS$ x S§)n 
plus some appropriate boundary conditions is dual to N = 4 d = 3 + 1 U (N) super- 
Yang-Mills. It means that, possibly, for an appropriate description of the known 
physical interactions through superstring theories, the strings have to exist in an 
Anti-DeSitter spacetime. The Anti-DeSitter spacetime has a negative cosmological 
constant, therefore it seems worthwhile to study spacetimes with a negative cosmo- 
logical constant if one wants to understand more about a 'possible' initial state of 
our primordial Universe. Of course, after that initial state the Universe would have 
to undergo a transition where the cosmological constant would change sign. Besides 
that, several important theoretical results and predictions in quantum cosmology 
have been obtained with a negative cosmological constant J23l24|25 l26l Considering 
a subset of all four-dimensional spacetimes with constant negative curvature and 
compact space-like hypersurfaces, S. Carlip et al showed how to compute the sum 
over topologies leading to the no-boundary wave-function.^EH These spacetimes 
are curved only due to the presence of a negative cosmological constant. In Ref. [53] 
it was shown how to obtain a vanishing cosmological constant as a prediction from 
the no-boundary wave-function and in Ref. [M]it was shown how to obtain predic- 
tions about the topology of the Universe from the no-boundary wave-function. We 
may also mention the result in Ref. 1251 where the WKB no-boundary wave-function 
of a homogeneous and isotropic Universe with a negative cosmological constant was 
computed. Due to the regularity condition imposed upon the spacetimes contribut- 
ing to the no-boundary wave- function, it was shown that only a well defined, discrete 
spectrum for the cosmological constant is possible. It was also found that among 
the spacetimes contributing to the wave function, there were two complex conjugate 
ones that showed a new type of signature change. 

The present paper is organized as follows. In Sec. [2j we introduce the classical 
model and obtain the appropriate Friedmann equation. With the aid of the poten- 
tial curve coming from this equation, we comment on the general behavior of the 
classical solutions. In Sec. [3j we quantize the model by solving the corresponding 
Wheeler-DeWitt equation. The wave-function depends on the scale factor a and on 
the canonical variable associated to the fluid, which in the Schutz variational for- 
malism plays the role of time (T). We separate the wave- function in two parts, one 
depending solely on the scale factor and the other depending only on the time. The 
solution in the time sector of the Wheeler-DeWitt equation is trivial, leading to an 
imaginary exponential of the type e~ lEr , where E is the system energy and r = — T. 
The scale factor sector of the Wheeler-DeWitt equation gives rise to an eigenvalue 
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equation. We find approximate solutions. In Sec. [4j we construct wave packets from 
the eigenfunctions and compute the time-dependent, expectation value of the scale 
factor. We find that the expectation value of the scale factor shows bounded oscilla- 
tions. Since the expectation value of the scale factor never vanishes, we confirm that 
this model is free from a big bang singularity, also, at the quantum level. Finally, in 
Sec. [3 we summarize the main points and results of our paper. 

2. The Classical Model 

The present Friedmann-Robertson- Walker cosmological model is characterized by 
the scale factor a(t) and has the following line element, 

ds 2 = -N 2 (t)dt 2 + a 2 (t) + rW) , (1) 

where df2 2 is the line element of the two-dimensional sphere with unitary radius, 

N(t) is the lapse function and we are using the natural unit system, where h = c = 

8ttG = 1. In this model the spatial sections are some closed three-dimensional solid 

,,[971 

with negative constant curvature, locally isometric to H J y^- Li The matter content 
of the model is represented by a perfect fluid with four- velocity C/ M = 5q in the 
comoving coordinate system used, plus a negative cosmological constant (A). The 
total energy-momentum tensor is given by, 

Tftu = (w + yjU^Uy - pg^ v - Ag^ , (2) 

As mentioned above, here, we assume that p = w, which is the equation of state for 
stiff matter. 

Einstein's equations for the metric ([TJ and the energy momentum tensor ([2]) are 
equivalent to the Hamilton's equations generated by the total Hamiltonian N(t)1-{., 
where, 



n = -^+3a + Aa 3 + ^. (3) 
12a a 6 



The variables p a and pr are the momenta canonically conjugated to the variables 
a and T, respectively. 

The classical dynamics is governed by the Hamilton's equations, derived from 
the total Hamiltonian iV(i)"H Eq. ([3]). In the gauge N(t) = a(t), they are, 

, = dNH _ _Pa 
dp a T> ' 



Pa= -dNK = -Qa-AK^ + 2^ 
dNU _ 1 



T = 



<5a~ a" 

(4) 



dpr ' 



„• _ dNU _ n 
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An Early Universe Model with Stiff Matter and a Cosmological Constant 5 

Where the dot means derivative with respect to the conformal time r = Nt, which 
in the present gauge is equal to at. 

We also have the constraint equation T~L = 0. It gives rise to the Friedmann 
equation, 



+ V c (a)=0, 



where the potential V c (a) is equal to, 

V c (a) = - 



Aa 4 
~3~ 



Pt 
3a 2 ■ 



(5) 



(0) 



For the present situation where A < 0, we have bounded solutions. The classical 
solutions are bouncing ones, free from an initial singularity, because we are con- 
sidering a stiff matter perfect fluid with negative energy density. Those results can 
be directly seen from the potential expression. A particular example of V c {a), for 
A = —0.1 and pt — —180, is given in Figured] 

3. The Quantum Model 

We wish to qua ntize the mo del following the Dirac formalism for quantizing con- 
strained systems .EEESEQEn] p- rs j. we introduce a wave-function which is a function 
of the canonical variables a and T, 



* = *(a,f). 



(7) 



Then, we impose the appropriate commutators between the operators a and T and 
their conjugate momenta p a and fir. Working in the Schrodinger picture, the oper- 
ators a and T are simply multiplication operators, while their conjugate momenta 



V(a) 




12345 6 789 10 



Fig. 1. V c (a) for A = 



—0.1 and pt = 



-180. 
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are represented by the differential operators 

Pa ~i Ta , PT -+ ■ (8) 

Finally, we demand that %, the superhamiltonian operator corresponding to ([3]), 
annihilate the wave-function 1 J> Eq. ([7]), which leads to Wheeler-DeWitt equation 

-l^£-( 3a2 + Aa4 ))* (a ' T)= V^^ < 9 > 

where the new variable r = —T h as b een introduced. 

The operator % is self-adjoint ^2 with respect to the internal product, 

da^^(a,T)*^(a,T), (10) 



cr 



if the wave- functions are restricted to the set of those satisfying either ^(0, r) = 
or ^'(0, r) = 0, where the prime / means the partial derivative with respect to a. 
The Wheeler-DeWitt equation ([9|) may be solved by writing ^(a,T) as, 

*{a,T)=e- iEr r,{a) (11) 

where 77(a) depends solely on a. Then 77(a) satisfies the eigenvalue equation 

-^^- + V(a) V (a) = UE±r)(a), (12) 
da a z 

where the potential V(a) is given by 

V{a) = -36a 2 - 12Aa 4 . (13) 

In the same way as in the classical regime, the potential V(a) Eq. ([T3]) gives 
rise to bound states. Therefore, the energies E, Eq. (fT2")l , of those states will form a 
discrete set of eigenvalues E n , where n — 1,2,3,.... For each eigenvalue E n there will 
be a corresponding eigenvector rj n (a). The general solution to the Wheeler-DeWitt 
equation © is a linear combination of all those eigenvectors, 

00 

^(a,T) = ^C„77 n (a)e- li5 " r , (14) 

n=l 

where C„ are free coefficients to be specified. 

We are going to use the Galerkin or spectral method (SM), in order to solve the 
eigenvalue equation (fl"2"|) . This met hod is well presented in Ref.[33]and it has already 
been used in quantum cosmology! *^ ^ In the SM, one must choose orthonormal 
basis functions in order to expand the solution to the eigenvalue equation. The 
solutions to the present eigenvalue equation (fT2|) must fall sufficiently fast for large 
scale factor values (a). It means that we must restrict the initial infinity domain of 
our variable a, to a finite domain. Say, < a < L, where L is a finite number that 
has to be fixed. Here, we shall consider solutions satisfying the condition ^(0, r) = 0. 
Putting together all the above properties of the solutions to Eq. (fT2|) , it is convenient 
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An Early Universe Model with Stiff Matter and a Cosmological Constant 7 

(but not mandatory) to choose our basis functions to be sine functions. Therefore, 
we may write r] n (a) in Eq. (|12[) as, 

n=l 

where the A^s will be determined by the SM. In the same a domain, we may also 
expand, in the same basis, the other two important functions of a appearing in Eq. 



V(a)j] n (a) = ^2B n J- sin (—J , (16) 

n=l 

- U„( a ) = ^B;J-sin — , (17) 



\ L I 

Where V(a) is given in Eq. (IT5|) and the coefficients B n and B' n can be easily 
determined. They are determined with the aid of Eq. (fT5")) and the fact that the 
basis functions are orthonormal. The coefficients B n and B' n are given by, 

oo 

(18) 

m— 1 

oo 

K = £ C 'm,nAm, (19) 



m— 1 



where, 



^.^i^sin^n^sin^rfa, (20) 



^/ 2 f L . /mira\ f 12\ . /mra\ , 

^,n = Z / o Sln (^)Uj Sln ( — ^ 

Introducing the results Eqs. (fT5)) - (l2"T1) in the eigenvalue equation (TT2"j) and using the 
fact that the basis functions are orthonormal, we obtain, 

n oo oo 



L J 



(22) 

m— 1 m— 1 

In order to derive some numerical results we must fix a finite number of basis 
functions, in other words, a finite number for the maximum value of the summation 
indices. Let us call this number N. The greater the value of N, the closer our results 
will be to the exact ones. We shall be restricted by our computational resources. 
Equation (|22|) . may be written in a compact notation as, 

D'^DA =EA, (23) 
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where D and D' are N xN square matrices and their elements are obtained from Eq. 
(|2"2"j) . The solution to Eq. (|2Uf gives the eigenvalues and corresponding eigenfunctions 
to the bound states of our quantum cosmology model. 

It is important to mention that the most correct form of 77(a) in the limit when 
a — > is not given by Eq. (I15[) . In order to obtain that expression for 17(a), one 
has to introduce the ansatz 77(a) — Ca a (where C is a constant and a is a number 
to be determined) in Eq. (|T2|) . After that, one has to discard the terms which 
have as coefficients the cosmological constant (A) and the curvature of the spatial 
sections (k = —1). The main motivation to discard those terms is that they are 
proportional to a Q+4 and a a+2 and should be less important, in the limit a — > 0, 
than the term whose coefficient is an algebraic equation for a. Finally, one solves 
the resulting equation imposing that the coefficient of the only remaining term, 
proportional to a a ~ 2 , vanishes. This gives rise to a second order algebraic equation 
for a, which solution satisfying the boundary condition 77(0) = is given by: a = 
0.5 + y/0.25 + 12|£'|. Therefore, the most correct form for 77(a) in the limit when 
a — > is given by, 

77(a) = Ca 0.5+V0.25+12|£|_ (24) 

We notice that it cannot represent the correct expression for 77(a) in the limit a — > 00 
because it diverges in that limit. That solution is already known in the literature!^ 

4. Energy Spectrum, Wave Packet and Mean Value. 

In this section we will solve Eq. (|12l) using the SM. In order to choose the number 
of basis functions N and the values of A, we performed the following numerical 
procedures. First of all, in order to choose the values of A, we solved numerically, 
the eigenvalue equation (|23[) for several different values of A and fixed values of 
N and L. We noticed that, although, we are free to choose any value of A, the 
results accuracy for small absolute values of A is better than for large absolute 
values of A, for a given number of basis functions N. This means that, if we use 
a large absolute value of A, we need to increase N to obtain the same accuracy 
of the case with a small absolute value of A. This, of course, would increase the 
computation time. Therefore, we have decided to use small absolute values of A. 
Taking in account this considerations, we choose A = —0.1 in the rest of our paper. 
In order to choose the value of A, we solved numerically, the eigenvalue equation 
(f23| for several different values of N and fixed values of A and L. Then, we compared 
the eigenvalues coming from the different choices of N. We noticed that only some 
of the first eigenvalues, for each different choice of N, remained the same up to a 
satisfactory accuracy. Therefore, we have decided to use N = 100 and take only 
the first eighteen eigenvalues to construct the wave packet. The accuracy of the 
eigenvalues, in this case, is in the tenth digit after the dot. Finally, it is important 
to men tion that based on comparisons with results of other models studied with the 
SMpSlwe chose L = 6, in the present case. 
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An Early Universe Model with Stiff Matter and a Cosmological Constant 9 

Now, using all these values: N = 100, A = —0.1 and L = 6, in the determinant 
constructed from Eq. (l23l) . we solve it and obtain the first 100 energy eigenvalues 
for the present case. From those, we take the first 18 energy levels and list them in 
Table [U 



Table 1. The eighteen lowest energy levels for a FRW model with k = — 1, A = —0.1, a 
stiff matter perfect fluid (p = p), N = 100 and L = 6. 



Ex = -380.2201284331828 


Ei = -342.1147751350869 


E 3 = -305.9147225014253 


E 4 = -271.6319016521779 


E 5 = -239.2791871064332 


E 6 = -208.8705235210961 


E 7 = -180.4210774184946 


E s = -153.9474207233139 


E 9 = -129.4677552918257 


E 10 = -107.0021912238143 


E u = -86.57309703051190 


E 12 = -68.20554802882278 


E13 = -51.92791264850499 


E 14 = -37.77263895532904 


E 15 = -25.77734426711615 


E 16 = -15.98638945709977 


E 17 = -8.453288958554614 


E 1S = -3.244733126937446 



In order to give an idea how the energy spectrum depends on A, we have con- 
structed the curve of the fundamental energy level E\ versus A. It is given in Figure 
EJ We notice that E\ decreases when one increases A. 

It is important to mention that even though the expression of r\(a) given by 
Eq. (|15j) is not the most correct one, in the limit a — > 0, we were able to determine 
numerically that it converges rapidly to zero in that limit. In particular, for the first 
energy level the corresponding eigenfunction oscillations about the zero value, in a 
region very close to it (0 < a < 0.2), are of the order of 10 -12 . For the eighteenth 
energy level the corresponding eigenfunction oscillations about the zero value, in 
a region very close to it (0 < a < 0.04), are of the order of 10 -9 . For the other 
eigenfunctions the oscillations about zero, in regions very close to zero, have values 
between the ones given above. 

Next, we construct the wave packet \l/(a,T) Eq. (TT41 . with the aid of rj n (a) Eq. 



El 01 



70- 
80- 
90- 



Fig. 2. Dependence of Ei with A, for N = 100 and L = 6. 
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(|T5)) and the energy levels in Table [TJ Our numerical study showed that, the energy 
eigenfunctions in the linear combination Eq. (|14[) are orthonormal and ^(a^r) Eq. 
(fl"4|) . has a constant norm. In the linear combination Eq. (TTTl) . we set C n equal to 
one, for the first eighteen values of n and C n equal to zero for the other values 
of n. In Figure [3j we show, as an example, the modulus squared of a wave packet 
constructed with the first eighteen energy levels, given in Table [TJ for r = 1000, 
L = 6 and A = -0.1. 

Finally, using the wave packet ^(a, r) we compute the mean value for the scale 
factor a, according to the following expression, 

Jo a 2 \^(a,T)\ 2 da 

We computed this quantity for many different time intervals. For all this different 
values, we observed that (a) performs bounded oscillations and never assume the 
zero value. Therefore, we confirm that these models are free from singularities, also, 
at the quantum level. As an example, we show in Fig. 2] the mean value computed 
with the wave packet constructed with the first eighteen energy levels, given in 
Table UJ for the interval from r = until r = 1000, L = 6 and A = -0.1. 



5. Conclusions. 

In the present paper, we quantized a Friedmann-Robertson- Walker model in the 
presence of a stiff matter perfect fluid and a negative cosmological constant. We used 
the variational formalism of Schutz. The model has spatial sections with negative 
constant curvature. The quantization of the model gave rise to a Wheeler-DeWitt 
equation, for the scale factor. We found the approximate eigenvalues and eigenfunc- 
tions of that equation by using the Galerkin or spectral method. In the present 




Fig. 3. Modulus squared of a wave packet constructed with the first eighteen energy levels for 
r = 1000, L = 6 and A = -0.1. 



November 4, 2011 0:44 WSPC/INSTRUCTION FILE 

G-Oliveira-Neto'corrected 



An Early Universe Model with Stiff Matter and a Cosmological Constant 11 

work, we considered only the negative energy eigenvalues and their corresponding 
eigenfunctions. This choice implies that the energy density of the perfect fluid is 
negative. A stiff matter perfect fluid with this property produces a model with a 
bouncing solution, at the classical level, free from an initial singularity. After that, 
we used the eigenfunctions in order to construct wave packets and evaluate the time 
dependent, expected value of the scale factor. We found that the expected value of 
the scale factor evolve with bounded oscillations. Since the expectation value of the 
scale factor never vanish, we confirm that this model is free singularities, also, at 
the quantum level. 
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